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Perturbative QCD results in the MS scheme can be dramatically improved by switching to a 
scheme that accounts for the dominant power law dependence on the factorization scale in the 
operator product expansion. We introduce the "MSR scheme" which achieves this in a Lorentz and 
gauge invariant way. The MSR scheme has a very simple relation to MS, and can be easily used to 
reanalyze MS results. Results in MSR depend on a cutoff parameter R, in addition to the /i of MS. 
R variations can be used to independently estimate i) the size of power corrections, and ii) higher 
order perturbative corrections (much like (x in MS). We give two examples at three-loop order, the 
ratio of mass splittings in the B*-B and D*-D systems, and the Ellis- Jaffe sum rule as a function 
of momentum transfer Q in deep inelastic scattering. Comparing to data, the perturbative MSR 
results work well even for Q ~ 1 GeV, and power corrections are reduced compared to MS. 



Introduction and Formalism 

The operator product expansion (OPE) is an impor- 
tant tool for QCD. In hard scattering processes two im- 
portant scales are Q, a large moment transfer or mass, 
and Aqcd, the scale of nonperturbative matrix elements. 
The Wilsonian OPE introduces a factorization scale A*, 
where Aqcd < < Q, and expands in Aqcd/Q- Con- 
sider a dimcnsionless observable a whose OPE is 



C w (Q,A/)C(A / ) + C , r(Q,A/) 



W) 



Q" 



(i) 



The Cq^ are dimensionless Wilson coefficients contain- 



ing contributions from momenta k > A^ with pertur- 
bative expansions in a s , and 6 Q V 1 = (0o,i)w are non- 
perturbative matrix element with mass dimensions and 
p, containing contributions from k < Af. If Cqi(Q,A/) 
are expanded they contain an infinite series of terms, 
(A f /Q) n , modulo ln m (A / /Q) terms, and this reflects the 
fact that C(5^ only include contributions from momenta 
k> N . The Wilsonian OPE provides a clean separation 
of momentum scales, but can be technically challenging 
to implement. In particular, it is difficult to define A^ and 
retain gauge symmetry and Lorentz invariance, and per- 
turbative computations beyond one-loop are atrocious. 

A popular alternative is the OPE with dimensional 
regularization and the MS scheme, which preserves the 
symmetries of QCD and provides powerful techniques for 
multiloop computations. In this case Eq. (Q) becomes 



Qp 



(2) 



where /z is the renormalization scale and bars are used 
for MS quantities. In MS the C\ are simple series in a s , 



C\(Q, M )=l + 5> n (g: 

n—1 



(4tt) 



(3) 



with coefficients b n (^/Q) = J2k=o ^ nk mfe (/VQ) contain- 
ing only In fi/Q. We will always rescale a and the matrix 



elements Qi such that Ci = 1 at tree level. In MS all 
power law dependence on Q is manifest and unique in 
each term of Eq. (0) . Also simple renormalization group 
equations in /i, like dhiCo{Q, /i)/rfln/z = j[a 8 (/j,)], can 
be used to sum large logs in Eq. (||) if Q ^> Aq CD . 

Cj v (Q,A^) and Cj(Q,fi) are related to each other in 
perturbation theory, so Eqs. (|l|) and @ are just the same 
OPE in two different schemes. The renormalization scale 
fi in MS plays the role of AA This is precisely true for log- 
arithmic contributions, In /i In A^, and here the Wilso- 
nian picture of scale separation in (7j and §i carries over. 
However, the same is not true for power law dependences 
on A^ . MS integrations are carried out over all momenta, 
so the Ci actually contain some contributions from arbi- 
trary small momenta, and the Qi have contributions from 
arbitrary large momenta. For the power law terms there 
is no explicit scale separation in MS, and correspondingly 
no powers of /i appear in Eq. (|[). While this simplifies 
higher order computations, it is known to lead to facto- 
rial growth in the perturbative coefficients. For Co, one 
has b n+1 (n/Q) ~ (ii/Q)Pn\ [2p /p] n Z at large n §, for 
constant Z. In practice this sometimes leads to poor con- 
vergence already at one or two loop order in QCD. This 
poor behavior is canceled by corresponding instabilities 
in 6\, and is referred to as an order-p infrared renor- 
malon in Co canceling against an ultraviolet renormalon 
in 6~i 1^, ||, U . The cancellation reflects the fact that the 
MS OPE does not strictly separate momentum scales. 

The OPE can be converted to a scheme that removes 
this poor behavior, but still retains the simple compu- 
tational features of MS. Consider defining a new U R- 
scheme" for Co by subtracting a perturbative series 

C (Q,R,n) = C Q (Q,n)-6C Q (Q,R,n), 



5C (Q,R,iJ,) = (^yj2 d 



(4tt) 



(4) 



with d n (fi/R) = J2k=o dnk ln fe (p/i?). If for large n the 
coefficients d n are chosen to have the same behavior as b n , 



2 



so d n+1 (n/Q) ~ (/^/i?)Pn![2/3 /p]"Z, then the factorial 
growth in Cq(Q,/x) and 5C${Q,R, ji) cancel, 



C (Q,R,fx) 



z. 



(5) 



Thus the i?-scheme introduces power law dependence 
on the cutoff, (R/Q) p , in Cq(Q, R, fi), which captures 
the dominant (A? /Q) p behavior of the Wilsonian . 
In practice this improves the convergence in Co even 
at low orders in the a s series. The dominant effect of 
this change is compensated by a scheme change to 9i, 
0iO) = 0i (i?, a«) - lQ p 5C (Q,R,n)}do{fi), and the new 
9± will exhibit improved stability. In the i?-scheme the 
OPE becomes 

a = C (Q, R, At)0o(M) + ^(Q, [i)9i(R, \i)IQ v 

+ C' 1 (Q^)9' 1 (R, f i)/Q p + ..., (6) 

where 9[ = [Q P 6C }9 and C[ = 1 - Ci ~ a s . Both 
Co and 9\ are free of order-p renormalons. The series in 
C{ 9[ is Borel summable. In all examples below 9o is also 
renormalon free. The above procedure may be repeated 
for higher renormalons and the higher power terms in 
the OPE indicated by ellipses, to improve the behavior 
of these terms as well. At the order at which we work, we 
will consistently set C\ = 1 and drop 9[ in the following. 

To setup an appropriate R-scheme it remains to de- 
fine the d n . In the renormalon literature such scheme 
changes are well known for masses. For OPE predictions 
a "renormalon subtraction" (RS) scheme has been imple- 
mented in Ref. ||. In RS an approximate result for the 
residue of the leading Borel renormalon pole is used to 
define the d n , which adds a source of uncertainty. 

For our analysis we define the "MSR" scheme for Co 
by simply taking the coefficients of the subtraction to be 
exactly the MS coefficients. In general it is more con- 
venient to use In Co rather than Co, since this simplifies 
renormalization group equations. Writing the series as 



mC (Q,/z) 



(47T) 



(7) 



with a n (p/Q) = J2k=o a ™ fc m /VQ we define the MSR 
scheme by the series 



In 



(47 



This definition still cancels the order-p renormalon for 
large n, as in Eq. (|^). It yields the very simple relation 



C (Q,R,fi) = Co (Q, At) [Co OR, AO] 



-(R/QY 



(9) 



which must be expanded order-by-order in a s (At) to re- 
move the renormalon. Thus the coefficient Cq(Q, R, /i) 
for the MSR scheme is obtained directly from the MS 
result. Note Co(Q, Q, p) = 1 to all orders. The appropri- 
ate p is obtained from the MS OPE by p=dimension(0i) — 



dimension(#o)- MSR preserves gauge invariance, Lorentz 
symmetry, and the simplicity of MS. 

The appropriate values for R in Eqs. (0,^,0) are con- 
strained by power counting and the structure of large 
logs in the OPE. The power counting 9\ ~ Aq CD im- 



plies 9\ ~ -^qcD' so f° r ^ ne matrix element we need 
R ~ A 4 ^ Aqcd (meaning a larger value where pertur- 
bation theory for the OPE still converges), which mini- 
mizes ln(Ai/AQco) and \n(/j,/R) terms in 9i(R, fi, Aqcd)- 
On the other hand, Cq(Q, R, At) has ln(At/Q) and \n(fi/R) 
terms, and for R ~ Aqcd no choice of At avoids large logs. 
For R ~ a 4 ~ Q we can minimize the logs in Cq(Q, R, A 1 ), 
but not in 9i(R, /i, Aqcd)- When the OPE is carried out 
in MS this problem is dealt with using a At-RGE to sum 
large logs between Q and Aqcd- For MSR we must use 
devolution, an RGE in the R variable J(|. The appro- 
priate R-RGE is formulated with fi = R to ensure there 
are no logs in the anomalous dimension. For Co, 

R-?LteC (Q,R,R)=j[a 3 (R)}- (|) P 7 [a s (i?)], (10) 



dR 

where ^y[a s ] 



E^oTnMWfTr]'^ 1 and 7 [a s ] = 
S^o 7n[ a s(-R)/47r]™ +1 are the MS and R anomalous di- 
mensions. Here 7n _i = pa n0 - 2 ^m=i m a m0 fin-m-i 
and we are using the MS /3-function fid/dfia s (/i) — 
-a^)/(2nY£^ =0 /3 n [a s ( f j,)/^} n - The choice in Eq. (§) 
keeps Eq. (|10|) simple. In cases where 7 is absent we 
expect Eq. (p_0[) to converge to lower scales due to the 
(R/Q) p factor multiplying 7. For Ri > R the solution 
ofEq. © is [U^U^RuRa)} 

Co{Q, Ro, Ro) = Co{Q, Rx,Rx) Ur(Q, Ri,Ro) , (11) 

where is a usual MS evolution factor and Ur is the 
i?-evolution. For p — 1 the complete solution for Ur was 
obtained in Ref. ||. It is straightforward to generalize 
this to any p. At N fc+1 LL order the (real) result is 



, k(k) k 

U R (Q,R 1 ,R ) = expl (^£) P £]s 



j i-npbx 



x p^lTi-ph-j^h) - Ti-ph-j^h)] \, (12) 



with r(c, t) the incomplete gamma function and £0,1 



-27r/(A,a a (i2o,i)). A 



.0) 



and A 



(2) 
QCD 



<>( D — ^ e '' ^"QCD 



Re\-t) b \ 



Re t {—t) bl e ^l 1 are evaluated at a large ref- 
erence R with t = -2tt/ (0 o a s (R)), and Si = 1 /{20l), 
h = {01 - /3 /3 2 )/(4/34), and b 3 = {(3\ - 2/3 /3i/3 2 + 
0o03)/@Po)- Defining 7n = 7 „/(2/3 ) n+1 the coefficients 
of Ur needed for the first three orders of i?-evolution are 

So = 70, Si = 71 - (6i+p6 2 )7o , (13) 

5 2 = 72 - (b 1 + P b 2 )ji + [(l+ph)b2+p(pbl+b 3 )/2]j . 

Eq. (|) becomes C (Q,R 1> R 1 )U R {Q,R 1 ,R )U li (R 1 ,R ) 
x# (-Ro) + 9i(Rq,Rq)/Q p , and this result sums logs be- 
tween R± r~j Q and Rq ~ Aqcd- This gives natural R 
scales for coefficients and matrix elements in the OPE. 
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Heavy Meson Mass Splittings in MSR 

The MS OPE for the mass-splitting of heavy mesons, 



A? 



' H 



'H 



for H 



B.D. is Am 2 H 



<?g(™q,m)m g M + EiCi{mQ,n)2fi(n)/{3mQ) + 
C(AQ CD /mg), where tuq = m b or m c . 



/'G 



Here 



G^h v \B v )/3 is the matrix ele- 
ment of the chromomagnetic operator, and p\ for 
i = ttG, A, LS, AG are 0(Aq CD ) matrix elements [Q, 
with p| g (m) — (3/2)A/1q(^). At the order of our 
analysis tree level values for the Cj suffice, so with 
£ p (/x) = (2/3) [ficfa) + fijji) - pl s (fi) + p| G ( M )] , 



Am 2 H = C G (m Q ,p) n G (p) + Y, p (n)/m Q 



(14) 



Taking the ratio of mass splittings r = Am 2 B / Am|, gives 



C G (m b ,/i) 

C G {m c ,pL) p 2 G (p)\m b m< 



M ( 



(15) 



The first term in this OPE gives a purely perturba- 
tive prediction for r. G G is known to suffer from an 
0(Aqcd/wq) infrared renormalon ambiguity |7j, with 
a corresponding ambiguity in The three-loop 

computation of Ref. || yields, r = 1 — 0.1113| Qs — 
0.0780| Q 2 — 0.0755| Q 3 at fixed order with fi — m c , and 

r = (0.8517) LL + (-0.0696)anll + (-0.0908)annll in 
RGE-improved perturbation theory, with no sign of con- 
vergence in either case. In MS these leading predictions 
are unstable due to the p = 1 renormalon in Co- 

Lets examine the analogous result in the MSR scheme 



Am 2 H = C G (m Q , R, m)Mg(m) 



m Q 



Since p = 1 the MSR definition in Eq. (H) gives 



CG{mQ,R,n) = C G (m Q , fi)[C G (R, fi)] 



— R/rriQ 



(16) 



(17) 



where C G (m, fj,) is obtained from Ref. B and we expand 
in a s (fi). The OPE in MSR at a scale i?o ^ Aqcd gives 



r = 



C G {m b , Rq, i?o) S p (iio, -Ro) / 1 1 
C G (m c ,i? ,-Ro) Hg( R o) ^m b m, 



(18) 



Large logs in Cg(toq, Rq, Rq) can be summed with the R- 
RGE in Eqs. (|lT]— |l3 ) . For simplicity we integrate out the 
b and c-quarks simultaneously at a scale R\ ~ yjm b ra c 3> 
i?o — Aqcd- This scale for i?i keeps ln(i?i/rnb )C ) small. 
With R-evolution and Ur from Eq. (O) we have 



(19) 



C G (m b ,Ri,R$U R (m b ,R 1 ,R ) 

T = 

C G (m c , Ri,R$U R (m c , R\,Rq) 

| Z p {R ,R ) / 1 1_\ 

H%(Ro) W m c /' 

This expression is independent of i?i and i?o- Order- by- 
order, varying i?i about ^Jm b m c yields an estimate of 
higher order perturbative uncertainties, much like vary- 
ing n in MS. For Rq the dependence cancels between the 
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FIG. 1: Perturbative predictions at leading order in 1/mq for 
the ratio r of the B-B* and D-D* mass splittings in the MSR- 
scheme (solid) versus MS (dashed). The Ro dependence of the 
solid red curve provides an estimate for the power correction, 
independent of the comparison with the experimental data. 
Neither Ri nor /i variation is shown in the figure. 

first term in r and the £ p power correction. In MSR the 
term E p (f?o,i?o) is ~ Aqcd ana - can ^ e positive or neg- 
ative. One may expect that there is a value of Ro where 
T. p (Rq, Rq) vanishes. Thus keeping only the first term 
in Eq. (03) and varying Rq > Aqcd yields an estimate 
for the size of this power correction. This technique goes 
beyond the dimensional analysis estimates used in MS. 

Fig. [j] gives perturbative predictions for r at differ- 
ent orders using the first terms in Eqs. ([l5 19) with 
m b = 4.7 GeV, m c = 1.6 GeV, a s {y/m b m c ) = 0.2627 
and the 4-loop /3-function. The solid lines are from the 
MSR scheme, plotted as functions of Rq. The dashed 
lines are the fixed order MS results with /i = ^/m b m c . 
The MSR results exhibit a dramatic improvement in 
convergence over MS for a wide range of Rq values. 
Varying Ri = ^m 6 m c /2 to 2^/m b m c at N 3 LL(MSR) 
gives Ar ~ ±0.008, which is a significant improve- 
ment over /i variation in the same range for N 3 LO(MS) 
where Ar ~ ±0.068. The MSR results converge to an 
Rq dependent curve, whose dependence cancels against 
^ p (Rq, Rq), so the residual R dependence provides a 
method to estimate the size of this power correction. The 
range Rq — 0.7 GeV to 1.2 GeV keeps Rq below m c and 
above Aqcd and yields 

r = 0.860 ± (0.065) Sp ± (0.008) pcrt . . (20) 

This estimate for the E p power correction in MSR is in 

good agreement with experiment, r Gxpt = 0.886 (D u d , 

B^l) and 0.854 (D { s * ] , B^). MSR achieves a conver- 
gent perturbative prediction for r at leading order in the 
OPE, and a 1/toq power correction of moderate size, 
~ 0.065, significantly smaller than the dimensional anal- 
ysis estimate of Aqcd(1/w c — l/m b ) ~ 0.2 in MS. 

Ellis-Jaffe sum rule in MSR 

In MS the Ellis-Jaffe sum rule p for the pro- 
ton in DIS with momentum transfer Q is M\(Q) = 
[C B (Q, t i)e B + 0)(Q,/i)ao/9] + 0i(/i)/Q 2 . Cb,o are 
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FIG. 2: Perturbative results for the Ellis- Jaffe sum rule in the 
MSR, RS, and MS schemes, at leading order in 1/Q. For all 
curves the one parameter, So, is fixed by data at Q ~ 5 GeV. 
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FIG. 3: Uncertainty estimates in the MSR scheme and MS 
scheme for the Ellis- Jaffe sum rule at leading order in 1/Q. 

known at 3 loops |Q. The two leading order terms are 
written so that both coefficients and matrix elements are 
separately ^-independent: 9b = gA/12 + as/36 is given 
by the axial couplings qa = 1-2694 and a§ = 0.572 for 
the nucleon and hyperon, while ao is a Q independent 
MS matrix element. d\ denotes all 1/Q 2 power correc- 
tions with their Wilson coefficients at tree level. The MS 
coefficients are affected by a p — 2 renormalon jllj , which 
is removed in the MSR scheme. Eq. (|J) gives [i = B, 0] 

d{Q, R, R) = Ci{Q, R)[Ci(R, R)}' r2 ^ 2 . (21) 



With devolution the MSR OPE prediction is 

Mi(Q) = [C B {Q,Ri,Rx)Ug{Q,R u Ro)e B (22) 
+ C (Q, R x ,Rx)U Q R {Q, R x , R )& /9] + 0i(flo, #o)/Q 2 , 

where U^'° are given by Eq. (|l^) with p = 2 and the cor- 
responding (a n o) B '° determine the appropriate (7 n ) B '°- 

Figures |],|] show perturbative predictions for the Ellis- 
Jaffe sum rule at leading power in 1/Q, compared with 
proton data from Ref. Q. We use a s U GeV) = 0.2282, 
and the 4-loop (3 with 4 flavors. In Fig. we show order- 
by-order results for the MS scheme at /i = Q, and for the 
resummed MSR scheme with R\ — Q and Rq — 0.9 GeV. 
We fix a = 0.141 so that MS and MSR agree with the 
data for Q ~ 5 GeV. MS agrees well with the data for 
large Q, but turns away at Q < 2 GeV and no longer con- 
verges. In contrast the MSR results still converge quickly 
and exhibit excellent agreement with the data over a wide 
range of Q values. The NLL MSR result already has the 
right curvature, and the NNLL and N 3 LL curves further 
improve the agreement. We also include predictions in 
the RS scheme with subtraction scale Vf = 1.0 GeV from 
Fig.3dof Ref. g, which improve slightly over the MS re- 
sults, but may not be capturing the dominant power law 
dependence on the factorization scale. In Fig. |^ we show 
uncertainties for three loop results in the MS and MSR 
schemes. The dashed red curve is the MS prediction, and 
the blue band estimates the higher-order perturbative 
uncertainties varying \x in the range /i mm (Q) < fi < 2Q. 
For Q > 1.5 GeV, [i mm = Q/2, while for Q < 1.5 GeV, 
(j™ in = 1.3Q/(1.1 + Q/(1 GeV)). The red solid line is the 
MSR prediction, the red band is the perturbative uncer- 
tainty from varying R\ in the same range as was done for 
/i in MS, and the green band estimates the 1/Q 2 power 
correction by varying Rq = 0.7 to 1.2 GeV. Fig. || implies 
-0.01 GeV 2 < 0i(i? o ,#o) ^ 0.01 GeV 2 in MSR, which is 
a much smaller power correction than the ~ 0.1 GeV 2 es- 
timate obtained from naive dimensional analysis in MS. 
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